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Abstract. The present paper aims at investigating the effect of rotation up to third 
order in the angular velocity of a star on p and g modes, based on the formalism 
developed by Soufi et al. (1998). Our ultimate goal is the study of oscillations 
of [3 Cephei stars which are often rapidly rotating stars. Our study shows that 
the third-order perturbation formalism presented by Soufi et al. (1998) must be 
corrected for some missing terms and some misprints in the equations. As a first 
step in our study of (3 Cephei stars, we quantify by numerical calculations the 
effect of rotation on the oscillation frequencies of a uniformly rotating zero-age 
main-sequence star with 12 Mq. For an equatorial velocity of 100 kms -1 , it is 
found that the second- and third-order corrections for (/, m) = (2, 2), for instance, 



are of order of 0.01% of the frequency for radial order n = 6 and reaches up to 
0.5% for n = 14. 
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1. Introduction 

Pulsating stars on the upper main sequence, and particularly S Scuti and (3 Cephei 
stars, are rapid rotators as well as being multimode pulsators. The ratio e = O/w of the 
rotation rate to the typical frequency of oscillations lo seen in these stars is no longer 
a small quantity as it is for e.g. the Sun. These stars have typically equatorial velocities 
~ 100 kms -1 , and oscillation periods from half to a few hours which implies e ~ 0.1, 
whereas for the Sun e ~ 10~ 4 . The effect of rotation on stellar structure and stellar 
oscillations is usually calculated through a perturbation analysis in which e is the small 
parameter. Whereas for the Sun a first-order perturbation analysis is sufficient, given the 
accuracy of observed oscillation frequencies, for these stars it is not. In order to achieve 
the full potential of asteroseismology for testing of models for upper main sequence stars 
a more careful treatment of the effect of rotation on oscillation frequencies is required. 

Rotation not only modifies the structure of the star but also changes the frequencies 
of normal modes. It removes mode degeneracy creating multiplets of modes. If the ro- 
tational angular velocity, VL, does not have any latitudinal dependence, and the rotation 
is sufficiently slow, the multiplets show a Zccman-likc equidistant structure. At faster 
rotation rates non-negligible quadratic effects in f2 cause the position of the centroid 
frequency of multiplets to shift with respect to that of a non-rotating model of the same 
star. 

Our long term goal is to study the oscillation properties of rapidly rotating (3 Cephei 
star. In this paper, we start with the study of rotation up to third order in the angular 
velocity of the star on p and g normal modes. To do this we use the third-order perturba- 
tion formalism according to Soufi et al. (1998), hereafter S98, but correct some misprints 
and missing terms in some of their equations. We carry out numerical calculations for 
the frequency corrections for a stellar zero-age main-sequence (ZAMS) model with mass 
M — 12 Mq, Mq being the solar mass. Section 2 gives a brief overview of previous 
relevant work on rotational perturbation theory for stellar oscillations. Section 3 shows 
the equations for the effect of rotation on oscillations up to third order and discusses 
the zeroth-order eigensystem. Section 4 presents corrections to eigenfrequencies and cou- 
pling coefficients. The numerical results are presented in section 5. Section 6 is devoted 
to concluding remarks. A summary of the differences between the present formulation 
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and that of b98 is given in Appendix A. Lengthy tormulae are collected m Appendix a. 
The hermiticity properties of an oscillating rotating system are discussed in Appendix C. 

2. Previous rotational perturbation analyses 

Simon (1969) and Saio (1981) studied the frequency corrections due to rotation up to 
second order for polytropes. Chlcbowski (1978) calculated the same corrections for white 
dwarf models. Gough & Thompson (1990) studied the effects of rotation and magnetic 
field on stellar oscillations up to second order. They investigated the linearized adiabatic 
oscillation equation in a rotating frame under the Cowling approximation and without 
viscous and resistive dissipative forces. They considered only an axisymmetric magnetic 
field, but allowed for its axis not to coincide with the rotation axis. They concluded 
that rotation and a magnetic field not only split the degenerate frequency multiplets but 
also shift the central frequency of each multiplet. The shift arises both from the direct 
effect of the perturbed inertial and Lorentz forces on the waves, and also because the 
unperturbed centrifugal and Lorentz forces change the structure of the star, which is 
no longer spherically symmetric. Gough & Thompson (1990) succeeded in formulating 
the differential equations governing the oscillatory motions in a form that is hermitian. 
The hermiticity of the oscillation equations for the rotating star was preserved by means 
of an appropriate mapping of each point in the distorted model to a corresponding 
point in the spherically symmetric stellar model. Although their main attention was on 
the effect of a magnetic field on the normal modes, they also performed some numerical 
calculations of the effect of rotation for three different, latitudinally independent, angular 
velocity profiles f2(r). They found that the effect of the second-order centrifugal distortion 
changes little with spherical harmonic degree, I, and can be approximated well by the 
asymptotic estimate. Also, for low I the second-order correction due to the advection 
term is negligible compared with the second-order centrifugal distortion and only for 
I > 50 are the two comparable in magnitude. 

Dzicmbowski & Goode (1992) derived a formalism for calculating the effect of differ- 
ential rotation on normal modes of rotating stars up to second order. They considered 
angular-velocity profiles with both radial and latitudinal dependency and found that 
at faster rotation non-negligible quadratic effects in 51 cause a departure from equidis- 
tant splitting. They also obtained generalized asymptotic formulae for g-mode splitting 
for which the Coriolis term is included in the zero-order treatment. These asymptotic 
results are relevant for white dwarfs and S Scuti stars. Dziembowski & Goode (1992) 
concluded that for solar oscillations the second-order effects arc dominated by distortion 
for I < 500. 

Soufi et al. (1998) extended the formalism of Dziembowski & Goode (1992) up to 
third order for a rotation profile that is a function of radius only. Their analysis has two 



advantages compared with previous investigations. By taking into account parts ot the 
effects of the Coriolis force in the zero-order system, the eigenvalue problem for stellar 
oscillations can be solved up to cubic order without having to solve successive equations 
for the cigenfunctions at each order. Also the usual m-degeneracy occurring in the absence 
of rotation is removed at the lowest order. S98 found that near-degenerate coupling due 
to rotation only occurs between modes with either the same degree I (and different radial 
orders) or with modes which differ in degree by 2. The first case involves modes in 
avoided crossings. The second case concerns modes that have close enough frequencies in 
the non-rotating model to be shifted into resonance if the rotation is sufficiently rapid. 
In general they showed that the total coupling comes from three distinct contributions: 
the Coriolis contribution, the non-spherically-symmetric distortion, and a coupling term 
which involves a combination of these two effects. 

Sobouti (1980) studied the normal modes of rotating fluids up to (3(f2 2 ). He argued 
that the p modes allow a perturbation expansion in fl, whereas this is not the case for 
the g modes. For, from a mathematical point of view, the condition for a perturbation 
series to converge is that the perturbing operator remain smaller than the unperturbed 
operator throughout the Hilbert space spanned by the normal modes; this condition is 
not met by g modes. 

Sobouti & Rezania (2001) considered the toroidal modes of rotating fluids, and showed 
that a) At 0(£l) the neutral toroidal motions of the non-rotating fluid organize them- 
selves into a sequence of modes with a definite (I, m) symmetry. Their radial degeneracy, 
however, persists at this order, b) Coupling of a given toroidal mode of (l,m) symmetry 
with (I ± 2, m) toroidal and with (I ± 1, m) poloidal modes, as well as the removal of the 
radial degeneracy come about at (9(fi 2 ). 

Result of calculation of frequency corrections up to third order were presented for 
models of (5-Scuti stars by Goupil et al. (2001), Goupil & Talon (2002), Pamyatnykh 
(2003), and Goupil et al. (2004). Daszynska-Daszkiewicz et al. (2002) studied the effects 
of mode coupling due to rotation on photometric parameters (amplitude and phase) of 
stellar pulsations. They reconfirmed the conclusion of S98 that the most important effect 
of rotation is coupling between close frequency modes of spherical harmonic degree, Z, 
differing by 2 and of the same azimuthal order, m. They presented some numerical results 
for a sequence of f3 Cephei star models with uniform rotation and for two- and three- 
mode couplings. Their calculations were carried out to cubic order in the ratio of rotation 
to pulsation frequency, according to the third-order formalism of S98. They concluded 
that due to the increasing effect of centrifugal distortion with the mode frequency, the 
coupling between acoustic modes is stronger than between gravity modes. 



6. I hird order perturbation formalism 

Following Unno et al. (1989), the equilibrium state of a rotating star can be characterized 
by a velocity field 

v = ft x r = f2rsin#e v , (1) 

where denotes the angular velocity. The rotation axis of the star lies in the 9 = axis 
of a spherical coordinate system (r, 9, <p). Q is assumed to be independent of latitude and 
can be written as: 

n = 0(r)e z = n[l + 77(r)]e z , (2) 

where Cl is the mean rotation rate. The stationary equation of motion in an inertial frame 
of reference is: 

-(v o -V)v o = ^+V0, (3) 
P 

where p, p, and 4> are the pressure, density and gravitational potential, respectively. One 
can show that the left hand side of Eq. (3) is equal to the centrifugal acceleration, F, in 
a corotating frame as follows: 

F = n x (fi x r) = rfl 2 sin6»e s = -(v • V)v , (4) 

where e z , and e s = sin0e r + cos9eg are unit vectors in cylindrical coordinates (s, ip, z). 
See Tassoul (2000). 

3.1. Equilibrium structure of rotating stars 

The stationary equation of motion, Eq. (3), is solved following Chandrasckhar (1933) 
and Chandrasekhar & Lebovitz (1962) by expanding the equilibrium quantities in terms 
of Legendre polynomials as 

f(r, 9) = f(r) + e 2 f 2 = f(r) + e 2 f 22 (r) P 2 (cos 9) , (5) 

where P 2 (cos9) = 3/2 cos 2 9 — 1/2 is the second Legendre polynomial and / can be p, 
p or 4>. For a rotation rate that is a function of r only, higher-order multipole moments 
(i.e., P4,Pe, ..) need not be considered. 

3.1.1. Spherically symmetric distortion 

The spherically symmetric part of the equilibrium structure can be obtained by substi- 
tuting Eq. (5) in the 9- independent part of the radial component of Eq. (3). The result 
is 



where the effective gravity is g c : 

g e = g - -rtf , (7) 



with 

„ d<j> GM r 



and M r is the mass confined within radius r. Eq. (6) is similar to the equation of hy- 
drostatic equilibrium of a non-rotating star but modified by including the spherically 
symmetric part of the centrifugal force. Note that the other equations governing the 
quantities of internal structure do not change; hence as in S98, following the approach 
described by Kippenhahn & Weigert (1994) a standard evolutionary code only modified 
according to Eq. (6) can be used to follow the evolution. Of course, in order to compute 
a model at a given age, the profile of the rotation rate is needed which requires the 
knowledge of its temporal evolution. 

3.2. Non-spherically symmetric distortion 

Substituting Eq. (5) in the ^-dependent parts of the radial and the tangential components 
of Eq. (3) and neglecting 0(e 4 ) terms yields the non-spherically-symmetric part of the 
equilibrium structure as: 

P22 = -pr 2 (Cl/e) 2 u 2 , (9) 
pr(n/e) 2 /dlnp 



P^ = f ^f^{^ + b 2 ) , (10) 
where 

u 2 =%n/e)- 2 + |(l + % ), 



(11) 

with T] 2 = i](ri + 2). The perturbed gravitational potential (f>22 satisfies the perturbed 
Poisson equation: 

< 12 > 

Equation (12) can be solved by numerical integration and with the appropriate boundary 
conditions: <\>2i oc r 2 at the centre of the star, and 022 oc r~ 3 at the surface of the star. 
The relations (9) to (12) are identical to Eqs. (77)-(80) (with k=l) of Dziembowski & 
Goode (1992) and also with Eqs. (15)-(17) of S98. 

3.3. Zeroth-order eigensystem 

Following S98, we include parts of the Coriolis and non-spherical distortion effects in 
the zero-order eigensystem. This yields cigcnfrcqucncics c^o of eigenmodes which are no 
longer m-dcgenerate, even at zero order. This way of building the zero-order eigensystem 



and the associated basis ot cigenmodes enables one to solve the eigenvalue problem up 
to cubic order without having to solve the successive equations for the eigenfunctions at 
each order. 

A zero-order mode £ is defined by 

£o = £ P o + e&i , (13) 

where the poloidal £ p o and toroidal £ t i eigenfunctions are characterized by a single spher- 
ical harmonic, as in the case of a non-rotating model, and are given by: 

C P o = r(yYr + ) 

o ( 14 ) 

= r ir ( Tl + lGr x v ^ r ™i + fl - lGr x v ^ y A) . 

where Vh = (§g e e + ^e^ p e 'p) 18 the horizontal part of the gradient operator and 
U)Q = ojq + 

According to Eq. (39) in S98 the eigenfrcquency uj n can be rewritten as: 

Uq = Wg 0) + U>1 , (15) 

where to^ is the usual normal-mode frequency of the system when rotation is absent, 
and 

wi = mQ{C h - 1 - Ji) . (16) 

The quantities Cl and Ji are given by Eq. (40) in S98 and references therein. Note that, 
following S98, we assume a dependence of the perturbations on time t and longitude ip 
as oc exp[i(m</? + cut)], such that prograde modes have m < 0. 

The frequency luq includes the first-order correction due to rotation, u>i, as well as 
the second-order correction due to spherically symmetric distortion. Also note that u>i is 
the usual first-order frequency shift due to the Coriolis force and can be rewritten as: 

^ = -jj n (0 [y 2 + ( A - l )* 2 - 2 H r4 P dr - ( 17 ) 

where / is the mode inertia, 

I = (C P o I ti P o) = f drpr^y 2 + Az 2 ) , (18) 
with A = 1(1 + 1). 

3.4. Poloidal eigenfunctions 

Following Unno et al. (1989), we define the dimensionless variables yt, v and w as 
g c r V P J g e r g c dr 



Then as in S98, the expression for the poloidal components is obtained as 

(20) 



r% = (V g - 3 + h 1 )y + ((- V g )y t + V g v 



A - 3y + Az , 



dr 



= (C r a 2 -A-^-jy+(A+l-U- X - fti)yt - Av , (21) 



dv 

r— = (l-U-x)v + w, (22) 

dw UA UV Q ( k UV q \ 

r-r- = z y+z — g -yt+ A-- + *)«;, 23 

dr 1 — <r r 1 — <r r \ 1 — ay / 

Az - Cj/t - /iiy = . (24) 



Here cr = R i u> /GM is the square of the dimensionless oscillation frequency and 
a^a + maa, a^-^L^, (25) 



1 dlnp dlnp _ 1 dlnp 

Txdlnr dlnr ' 9 " ^ dlnr ' 



rr dlnM r 2 ^ 2 ._ 

_ JrJRf^ C 

6 = jmjm] ' Cr = r^; ' (28) 

A = y g (y-y t +w) , a = 2m°^, (29) 

a 

2rtt 2 / n d\nn 2 \ 



A A Aa 

C = i 9 , hi = . (31) 

A - a C r a A - a 



Here M and R are the stellar mass and radius respectively and G is the gravitational 
constant. 

For radial oscillations (I = A = 0) the Poisson equation for £ p0 , Eq. (12), reduces to 

Thus <fr' can be eliminated from the problem; instead of Eq. (19) we define the dimen- 
sionless variable y t in this case as 

yt = —-, (33) 

9crp 
obtaining 

r df = iVg 3)y Vgyt ' (34) 



r^=[C r a 2 -A+ Y —jy + (A + l-U- X )yf (35) 
For radial oscillations there is no transverse component, i.e., z = 0. 



j. b. toroidal eigenrunctions 

Following again S98 the components r, f of the toroidal part, £ t i, for a mode k = 
(nfe,Zfe,mfe), are obtained as follows: 

T k+1 = T n fc ,Z fc + l,m fc 

^ (m. h3m fe ^ 



= i IK Afc+1 _ x 2P fe + 3m fc 4^d 



'(A h+1 -a k )\ ' l; 

T~k-1 = Tn k ,l k -l,m k 

= Va Pk „ \ hPk + 3m k g-d k ) , 

with 

ft = (l + r ] )(lk + 2){-yk+lkZk) , (38) 
A = (l + v)[h - l){y k + (Ik + l)z k ] , (39) 



A = V%Fr- (40) 



/3fe+i being defined similarly, but with l k replaced by l k + 1, and 



rffc = Wfe + Dk - Ca k z k j y J^: u 2 + hj + A fe u 2 • (41) 
4. Correction to eigenfrequencies and coupling coefficients 

As discussed in S98, in general the correction to the eigenfrequencies must take into 
account coupling between nearby modes. Thus we need to evaluate coupling coefficients 
between different modes. However, in the case of an isolated mode, k, the correction w c , 
to the eigenfrequency u = uj + w c , is obtained from 

Hkk _ nj 
W c — TTvT = itkk 



(0) 


= u, T +u, D +u; C + 0(e 4 ) , (42) 
where Hk q is the matrix of the coupling coefficients and is defined as 
H kq = £f± = r kq + V kq + C kq + 0(e 4 ) , (43) 

£Jkq 

and Jk q = \J uj^uj^q Iklq- Note that for an isolated mode each contribution lo t = T k k, 
up = V k k and ui c = Ckk to the frequency correction arises from the corresponding 
diagonal term of the interaction matrix H kk . 

In order to obtain expressions for the different contributions to the total coupling 
coefficients Hkq, see Eqs. (B1)-(B2) in S98, it is convenient to define a dimensionless 
radius x and density p: 

x = r/R , 

(44) 

P = P/Pc , 



2Iu) 



where p c is the value ot the density at the centre ot the (deformed) star in equilibrium. 
The second-order perturbation terms of pressure, density and gravitational potential are 
also made dimensionless by defining the following: 

P22 =P22/Pc , 

P22 = P22/Pc , ( 45 ) 
022 = 022/^c , 

where p c = p c (j) c and <ft c = R 2 Cl 2 . Oscillation frequencies are made dimensionless by 
making use of the dynamical time scale of the star: 



(0) 

(0) 

Q _ 



Of2 



^GM/R 3 

n 



(46) 



^GM/R 3 ' 
Dimensionless mode inertia terms are: 



J] 



Jk 



kq ~ Pc R*JmfW ^ 

I (0) (0)1 f 
= y^k^Oq hlq ■ 

4.1. Coriolis Contribution: T kq 

As in S98 the elements T kq are obtained as 

T t = ^2 

kq JGM/W (48 ) 

W X T^ 1 ) _L. X ^ _L X ^ 2 > 

= 1 qk = 0l k l q l +0l k l q +2+ kq + d lkl q -2-L qk , 

where * denotes complex conjugate; here the diagonal and off-diagonal terms are: 



2J k . 



dxpx x 



X (A k+1 T^ +1 T q+1 + K k _ 1 T^_ 1 T q - 1 ) 



Tkq = 2jtj\J Ax ~ 9X Afc - lf fe-i T 9+ 1 



Amcr, 



n 



<?Q k + O09 



^ dxpx 4 (l +j?)ffe_ 1 r 9+ i| . 



(50) 



We recall that some effects of distortion are already included in the toroidal zcroth- 
order system. Therefore distortion contributes indirectly through the toroidal part of the 
eigenfunctions, i.e., through the components r and f. 

Note that the case k 7^ q represents a coupling between near-degenerate modes. 
Equation (48) shows that the coupling occurs between modes either with same degree I 



(but generally dittcrcnt radial orders) or with degrees that diher by ±2. Ihc frequency 
correction, uj t , for a single mode can be obtained from the diagonal element Tkk as 



7r T T 

■l kk = c 2 + a 3 , 



y/GM/R? 

with the second-order contribution 
^ = \{J )(g) 2 j ( 1 d^x 4 (A fc+1 |r fc+1 



(51) 



(52) 



+ A fc _i | fk-i 



and the third-order contribution 
^= m(ai)(a)3/ CT.-l-Ji 



dxpx 4 x 



x (A fe+ i | T fc+ i | 2 +A fe _i | r k -i | 2 ) 

- / dar^ 4 (l + ?7)(| r fe+1 | 2 + | ffe_i | 2 ) 
Jo 



(53) 



Equations (48) to (53) are identical to Eqs. (B3)-(B7) of S98. 

Note that to separate the frequency correction into the second and third order in 
terms of e = ct^/ctq, one also needs to use the following approximations: 



Ml - ~ 1 _i_ I i _ 



i , m ~ 



l + m^(C L -l- Ji) 



(0) 



(54) 



^ * (f ) (Sf) 

cx 0(^ 2 ) + 0(ft 3 ) 



(«) 



(55) 



in which it is used that 

_ (0) /QN 

^*_?0- = m(C L - 1 - Jx) oc O I = ] = 0(1) . 
°n V". 



(56) 



4.2. Non-spherically symmetric distortion: T>k q 

Following again S98, the coefficients T>kq can be obtained in an explicitly symmetric form 
as follows: 

Vk « =v i k = 2jrJ d3x { {VlP2) v ■ c p° k v ■ ?p09+ 

where 

1 dtt 2 



(57) 



,2 dn 2 

P2W 1 + (l-P 2 )-r— 
3 dr 



W = y k y q Y h *Y q 

x (y k z q Y k *V H Y q + y q z k Y q V n Y k ) ■ V Hj P 2 , 



-r x 

3 dr 



(58) 



where Y k = Y lk ; here F 2 is the second Legendre polynomial and w\ is dehned by 

d /1\ dp 22 . Id ( p 22 \ dp 
ar \p J ar par \ p J or 
The relations (57) to (59) are identical to relations (B8) to (BIO) in S98. Note that V kq 
is non-zero only when m k — m q . 

The angular parts of the integral Eq. (57) can be evaluated analytically. Finally the 
result can be written in dimensionless form as follows: 

V t = . Vkq 
q \/GM/R 3 

= 6 lklg D + S lklq+2 (lp k f3 q+1 )D kq (60) 
+ Si k i q - 2 (^f3k+if3 q )D kq , 

with 

a 2 - f 1 

D = Q kk2 D kq + -J^- / dxpx 4 y k y q b 2 , (61) 

■Jkq JO 

and 

= (62) 

where the expressions for the different terms appearing in D kq , are given in Appendix B. 
The diagonal elements V kk reduce to the frequency correction lo° as 

'"'Tmm'^"^^- (63) 

where the second-order contribution is 
Oji\ 2 ( Jk 

and the third-order contribution is 



•f -(¥)(*) (^) t c 



u 

3 



(65) 

= m(C L - 1 - jo (a ) 

in which Z) is given by Eq. (61). 

Note that the relations (60) to (65) differ from the relations (B12) to (B21) in S98. 
The differences come from the fact that S98 have further performed integrations by 
parts in order to remove dln/5/dlnr = — (A + V g ) which sometimes can cause numerical 
oscillations of the integrand of Eq. (57). For the models which we have considered (see 
Sect. 5) no such numerical problem exists and we choose to compute the coefficient 
directly from the definition in Eq. (57). We furthermore note that the partial integration 
gives rise to surface terms which are neglected by S98; in fact, for models without an 
extended atmosphere they become considerable. We return to this subject in Sect. A.l. 

Also in the third-order terms of the relations, the approximation z ~ yt/Ca 2 ,, which 
is identical with the non-rotating case, instead of its exact relation Eq. (24) was used by 
S98. Indeed using the relation Eq. (24) adds higher-order corrections which in principle 
ought to be negligible. However, we find in our numerical results that it is not true and 
Eq. (24) must be used (see Sect. A. 2). 



4.5. Distortion and Loriolis coupling: L kq 

As in S98, the coupling coefficients C kq can be obtained as: 

= Si k i q C + 5 lklg+2 (lp k (3 q+1 )C kq ( 66 ) 

+ &l k l q -l{ 2 fik+\Pq)Ckq , 

where 

C = oT-^fi^ofc + ^o,) I' dx~px 4 C(l + V )x 

ZJ k q JO 

x[b 2 -(A + V g )u 2 }x ( 6? ) 
x [z k z q + Q kq2 (ykz q + y q z k + 3z k z q )] , 

Ckq= oy-^(o-ofc + ^o,) / da;^ 4 C7(l + ?7)x 

ZJ fc<? Jo (68) 

x[6 2 - (A + y g )M 2 ][yfcZ 9 + + 3z fc z 9 ] ; 

an expression for Q kq2 = J sin9d8d(f)Y k *Y q P 2 is given by Eq. (B18) of S98. The diagonal 
coefficients correspond to the contribution to the frequency correction ui c = C kk . This 
contribution has only terms of third order as follows: 

= j<?hf dx~px 4 C(l + V )[b 2 -(A + V g )u 2 ]x ( 69 ) 
x [z 2 + Qfc fc2 (22/z + 3z 2 )] . 
Equations (66) to (68) are identical to the relation (B22) in S98. However, the relation 
(69) differs from Eq. (B25) in S98 where again a partial integration has been carried out 
to eliminate the density derivative. 

Note that Eqs. (53), (65) and (69) show that for the case m — there are no third- 
order contributions to the frequency corrections. 

The expressions for the separate terms in Ji kq (see Eq. (43)), i.e., Eqs. (48), (57), and 
(66), are all explicitly symmetric and real- valued. Therefore one can conclude that the 
total coupling coefficient is symmetric, real-valued and therefore also hermitian, TL kq = 
Ti-qk- 

4.4. Close frequencies and mode coupling 

The existence of close frequencies can lead to large values of the perturbation terms to 
the extent that it invalidates the perturbative approach outlined in the previous sections. 
For such cases one needs to use a degenerate perturbation formalism. Following S98, the 
zeroth-order eigenfunction of an individual mode is written as a superposition of the 
degenerate eigenfunctions: 



where £ /c 1S the degenerate eigcnlunction ot the zcroth-ordcr system and the amplitudes 
Ak are expansion coefficients which can be normalized to 1. As same in S98, the solution 
of the eigenvalue problem for the eigenfunction £, which corresponds to an eigenvalue uj, 
satisfies the following linear system: 

AMfcfc + X] AqUkq = , (71) 
q^k 

where the components fikq are defined by 

[i k k = K - oj)2uj { ^I k , ^ 
Hkq = 2J kq Hk q , for q ^ k . 

One notes that \ik q — H q k because of the symmetric properties of Hkq- The eigenvalues, 
u, can be obtained from the existence condition of nontrivial solutions of Eq. (71), i.e., 

det[/iy] = 0, i,j = l,...,JV, (73) 

where det is the determinant and N is the total number of the near-degenerate modes. 
Here we consider the coupling between N — 2 modes. Then Eq. (73) reduces to 

(wi -w)(w 2 -w) - W? 2 =0 . (74) 

In this case the eigenfrequencies, lo — uj±, can be obtained trivially from the zeroth-order 
near-degenerate eigenfrequencies and the total corrections due to rotation and coupling 
are 



" ± =(^W +n 2 + 0(e*). (75) 



2 

2 y ~ v V 2 ~" 

Also, the normalized amplitude A = (Ai,A2) appearing in Eq. (70) can be derived from 
Eq. (71) as follows 

4 ±] =[i + (a4!W]- 1/2 

= [l + (Mi2//4?) 2 ]- 1/2 m , 

4 ±} - -(^Vm 12 )4 ±) 

= -(M12/M22 =) )4 ±) , 

where the quantities labeled by (±) are related to the eigenvalues u±. 

The coupling coefficients Hkq appearing in Eq. (72) show that the near-degenerate 
coupling only occurs between modes with the same m and with either the same degree I 
(and different radial orders) or with degrees which differ by ±2. 

5. Oscillations of a rapidly rotating B star 

In order to calculate the effect of rotation on normal modes, we consider a uniformly 
rotating, 12 M©, ZAMS model generated by the evolution code of Christensen-Dalsgaard 
(1982) (see also Christensen-Dalsgaard & Thompson 1999). The parameters of the model 
are listed at Table 1. 



Table 1. Stellar parameters of a rotating zero-age mam-sequence star m solar units. M, M CO nv, 
i?, i?conv 5 p c , and p c are the total mass, the mass of convective core, the radius, the radius of 
convective core, the central pressure and density, and denotes solar values; ctq and e are the 
dimensionless mean angular velocity and the perturbational expansion coefficient; Td yn , Trot, and 
Vrot are the dynamical time scale (free fall time), the equatorial period and velocity, respectively. 
For comparison note that T©d y n = 0.5 h, r Qro t = 25 d, V© ro t = 2 kms -1 . 

M — 12 M Q 
R = 4.32 R e 
p c = 3.72 x KT 1 p Qc 
= 1.38 x 1CT 1 
T dyn = y / R 3 /GM = 1.15 h 
Vrot = R& = R(Jn/T dyu = 100 kms- 1 



The behavior of some of equilibrium quantities of the model against fractional radius, 
x = r/R arc shown in Fig. 1; It shows that: 1) From the centre to the radius x = 
0.25, the star is in a convective regime where N 2 < and, outside of this radius is 
in a radiative regime where iV 2 > 0; 2) The spherically symmetric density p decreases 
smoothly from its maximum value to nearly zero at the surface. 3) The non-spherically- 
symmetric correction to the density P22 decreases to a negative maximum value at x = 0.3 
and then slowly increases to zero at x ~ 0.9; 4) The non-spherically-symmetric correction 
to the gravitational potential 4>2i is obtained as the numerical solution of the Poisson 
relation, Eq. (12), by a Runge-Kutta method with an adaptive step-size control. The 
absolute value of </>22 increases smoothly to its maximum value at the surface. 

5.1. Eigenfunctions 

The zero-order eigenfunctions are computed from the zero order eigenvalue problem with 
the pulsation code of Christensen-Dalsgaard (see Christensen-Dalsgaard & Berthomieu 
1991), modified according to Eqs. (20) to (24). 

In Fig. 2, the radial (y) and horizontal (z) components of the zero-order poloidal 
eigenfunctions as well as rp 1 ^ 2 £ ir /(R 2 pl^ 2 ) related to the radial energy density, where 
£ r = ry is the radial displacement, are plotted against the fractional radius x = r/R for 
the selected f and p modes with (l,m) = (2,2) and n — (0, 1,8). Note that in the case 
of radial oscillation (I = 0), the results are derived from the reduced set of differential 
equations composed of Eqs. (34) and (35). 

Figure 2 shows that: 1) The number of nodes coincides with the radial order n. 

2) At lower radial order, n, the energy is distributed throughout most of the star. At 
higher radial order, n, the displacement is large only in the outer part of the model. 

3) The magnitude of the horizontal component z decreases with increasing n (see also 



M conv = 0.34 M 

Rconv — 0.25 R 

p c = 4.83 x 10~ 2 p Qc 
: = Q/uj = Ofi/27r = 2.2 x 10~ 2 
T Iot = 2-KR/V Iot = 2.15 d 



bobouti (1980)). In other words, the higher-order p modes osculate radially rather than 
horizontally because they are acoustic and longitudinal waves. 

5.2. Eigen frequencies and corrections 

The zero-order cigcnfrcqucncy, erg, is derived from numerical solutions of Eqs. (20) to (24) 
by the modified pulsation code; note that using the eigensystem in Eqs. (20) to (24) the 
first-order frequency correction, u\, is implicitly included in (see Eqs. (15) and (16)). 
The second- and third-order Coriolis contributions, (crj, crj), the second- and third- 
order non-spherically-symmetric distortions, (erf, erf), and the third-order distortion 
and Coriolis coupling, crj, are derived from numerical integrations of Eqs. (52), (53), 
(64), (65), and (69). 

In Tables 2 to 5 the results of different contributions of frequency corrections due to 
effect of rotation up to third order are tabulated. In each table the selected p modes with 
(I, to) = (0, 0), (1, 1), (2, 1) and (2, 2) with n = (0, . . . , 14), (1, . . . , 15), and (0, . . . , 14) 
respectively, are considered. The fundamental modes are known with n = and labelled 
by f. The modes with n > 1 are labelled by pi, . . . , p„. 

Tables 2 to 5 show that: 1) The values of zero-order eigenfrequency, <Jo, increase when 
the radial order n increases. 2) For the case m = 0, the odd-order frequency corrections, 
i.e., in the present case those of first and third order, vanish. In this case only the 
second-order frequency corrections exist (see Eqs. (17), (52), (53), (64), (65), and (69)). 
3) For radial oscillation (I = 0) only the second-order Coriolis contribution exists and the 
second-order non-spherically-symmetric distortion which results from centrifugal force 
vanishes for a uniform rotation (see Eq. (64)). 4) In non-radial cases, at lower radial 
orders, the Coriolis contribution and non-spherically-symmetric distortion have the same 
order of magnitude. But at higher orders, the magnitude of the non-spherically-symmetric 
distortion becomes greater than Coriolis contribution. 5) With increasing n, the frequency 
correction due to distortion and Coriolis coupling increases and decreases alternatively. 

In Figs. 3 to 6 the results in Tables 2 to 5 are plotted. The plots show that for 
high-order p modes: 1) <to increases as a linear function of n, which is very similar to the 
pattern of vibration in a simple string; this also follows from asymptotic analyses of high- 
order acoustic modes (e.g., Vandakurov 1967; Tassoul 1980). 2) (crj, <rf) and (crj , cJp) 
have almost regular asymptotic behaviors as l/cr ocl/n and 1/<Tq oc 1/n 2 , respectively. 
See Eqs. (52), (64) and (53), (65). 4) The distortion and Coriolis coupling, crj, has no 
clear asymptotic relation (see Eq. (69)). 5) The total frequency shift corrected up to 
third order, cr c , for radial and non-radial modes decreases and increases, respectively, 
with increasing n; in the non-radial case this is caused by the dominance of the term erf 
which vanishes in the radial case. Note that for the case of (I, to) = (2, 2), in the diagrams 



oi a 2 and a break between t and pi appears. Ihis happens because the inertia tor 
the f mode is two order of magnitude greater than for pi. 

In Table 6, the results on third-order frequency corrections for the case of two near- 
degenerate modes, derived from Eqs. (75)-(76), are tabulated. As discussed above the 
coupling exists only for the two near-degenerate poloidal modes belonging to the same 
m but with I differing by ±2. 
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Fig. 1. Dimensionless equilibrium quantities including squared buoyancy frequency N 2 = 
(g e /r)A with A — (l/Ti)(dhip/dm;E) — (dlnp/d In x), spherically and non-spherically- 
symmetric contributions p and p22 to the density, and non-spherically-symmetric gravitational- 
potential contribution 022, in units of GM/R S , p c and R 2 Q 2 , respectively, against fractional 
radius x = r/R for a zero-age main-sequence star with M = 12M . 
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Fig. 2. Zero-order radial y (left) and horizontal z components (middle) as well as 
r P 1 ^ 2 £r/(R 2 Pc 7 * 2 ), related to the energy density (right), against fractional radius x = r/R for 
selected f and p modes with (I, m) = (2, 2) and n = (0, 1, 8) for the model described in Table 1. 
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Table 2. Values ol eigenlrequency cto (including the (J(ll) contribution Irom rotation; see 
Sect. 3.3), second- and third-order Coriolis contributions aj and aj , second- and third-order 
non-spherically-symmetric distortions a® and erf, third-order distortion and Coriolis cou- 
pling CT3 , and total frequency cr to t = <?o + cr c corrected up to third order, for p modes in 
a zero-age main-sequence star with M = 12M , (l,m) — (1,1) and n = (1, . . . , 15). Here 
<t c = aj + aj + 02 + o 3 + C3 is the total third-order frequency correction. All frequencies are 
in units of ^GM/R 3 = 2.42 x 1(T 4 s" 1 . 



Mode 


n 


o~o 


„T 
02 




-2° 






,T T 

^3 






-3° 






-3° 




ftot 


Pi 


1 


3.3700 


2.9403x10" 


-3 


-1.8726 


X 


10" 


-3 


-1.4269 x 


10" 


-4 


6.5305 x 


10" 


-5 


6.9849x10" 


-5 


3.3711 


P2 


2 


4.5499 


2.1727x10" 


-3 


-2.2274 


X 


10" 


-3 


-7.9057 x 


10" 


-5 


5.8518 x 


10" 


■5 


7.4479x10" 


■5 


4.5499 


P3 


3 


5.6388 


1.7589x10" 


-3 


-2.6952 


X 


10" 


-3 


-5.2684 x 


10" 


-5 


5.8734 x 


10" 


■5 


9.5457x10" 


■5 


5.6380 


P4 


4 


6.6669 


1.4833x10" 


-3 


-2.7774 


X 


10" 


-3 


-3.8095 x 


10" 


-5 


5.2159 x 


10" 


■5 


1.0610x10" 


-4 


6.6657 


P5 


5 


7.7061 


1.2849x10" 


-3 


-2.7012 


X 


10" 


-3 


-2.8820 x 


10" 


-5 


4.4457 x 


10" 


■5 


1.0074x10" 


-4 


7.7048 


P6 


6 


8.8212 


1.1286x10" 


-3 


-2.6485 


X 


10" 


-3 


-2.2305 x 


10" 


-5 


3.8528 x 


10" 


■5 


9.1212x10" 


■5 


8.8198 


P7 


7 


9.9704 


1.0030x10" 


-3 


-2.5225 


X 


10" 


-3 


-1.7663 x 


10" 


-5 


3.2778 x 


10" 


-5 


8.3252x10" 


■5 


9.9690 


P8 


8 


11.152 


8.9863x10" 


-4 


-2.0507 


X 


10" 


-3 


-1.4235 x 


10" 


-5 


2.4023 x 


10" 


-5 


7.8685x10" 


-5 


11.151 


P9 


9 


12.336 


8.0980x10" 


-4 


-5.5136 


X 


10" 


-4 


-1.1657 x 


10" 


-5 


5.8801 x 


10" 


■6 


7.8932x10" 


-5 


12.336 


P10 


10 


13.503 


7.3232x10" 


-4 


3.3305 


X 


10" 


-3 


-9.6733 x 


10" 


-6 


-3.2642 x 


10" 




8.4408x10" 


-5 


13.507 


Pll 


11 


14.635 


6.6613x10" 


-4 


1.1346 


X 


10" 


-2 


-8.1472 x 


10" 


-6 


-1.0309 x 


10" 


-4 


9.1742x10" 


-5 


14.647 


Pl2 


12 


15.741 


6.1527x10" 


-4 


2.3286 


X 


10" 


-2 


-7.0149 x 


10" 


-6 


-1.9741 x 


10" 


-4 


9.2581x10" 


■5 


15.765 


Pl3 


13 


16.851 


5.7802x10" 


-4 


3.7709 


X 


10" 


-2 


-6.1686 x 


10" 


-6 


-2.9945 x 


10" 


-4 


8.5150x10" 


■5 


16.889 


Pl4 


14 


17.971 


5.4718x10" 


-4 


5.8331 


X 


10" 


-2 


-5.4850 x 


10" 


-6 


-4.3534 x 


10" 


-4 


7.5827x10" 


■5 


18.029 


Pl5 


15 


19.065 


5.1966x10" 


-4 


1.1101 


X 


10" 


-1 


-4.9178 x 


10" 


-6 


-7.8261 x 


10" 


-4 


6.8432x10" 


-5 


19.176 



Table 3. Same as Table 2, lor p modes with ((, m) = (0, 0) and n = (1, . . . , 15). There is no <7 3 , 
<r!p, and a§ = for m = (see Eqs. (53), (65), and (69)). Also vanishes since I = m = and 
the model is a uniformly rotating star (see Eq. (64)). 



Mode 


n 


CT() 


„ T 

&2 




Ctot 


Pi 


1 


2.9941 


8.4876x10" 


-3 


3.0026 


P2 


2 


3.9460 


6.4400x10" 


-3 


3.9524 


P3 


3 


5.1853 


4.9008x10" 


-3 


5.1902 


P4 


4 


6.2697 


4.0532x10" 


-3 


6.2738 


P5 


5 


7.3206 


3.4714x10" 


-3 


7.3241 


P6 


6 


8.4203 


3.0180x10" 


-3 


8.4233 


P7 


7 


9.5645 


2.6570x10" 


-3 


9.5671 


P8 


8 


10.742 


2.3656x10" 


-3 


10.745 


P9 


9 


11.932 


2.1299x10" 


-3 


11.934 


PlO 


10 


13.112 


1.9382x10" 


-3 


13.114 


Pll 


11 


14.261 


1.7820x10" 


-3 


14.262 


Pl2 


12 


15.374 


1.6529x10" 


-3 


15.376 


Pl3 


13 


16.478 


1.5422x10" 


-3 


16.479 


Pl4 


14 


17.591 


1.4446x10" 


-3 


17.592 


Pl5 


15 


18.702 


1.3588x10" 


-3 


18.704 



Table 4. Same as Table 2, for f and p modes with (I, m) = (2, 1) and n = (0, . . . , 14). 



Mode 


n 


CT() 


„ T 

°2 










^3 






-3° 






-3° 






ftot 


f 





2.8330 


4.1956x10" 


-3 


7.3274 x 


10" 


-4 


-2.8916 


x 


10" 


-4 


-2.2495 


x 


10" 


■5 


-4.7041 x 


10" 


-5 


2.8376 


Pi 


1 


3.7712 


3.6956x10" 


-3 


1.3637 x 


10" 


-3 


-1.8491 


x 


10" 


-4 


-4.1312 


x 


10" 


■5 


-5.9259 x 


10" 


-5 


3.7760 


P2 


2 


4.9925 


2.9472x10" 


-3 


1.7496 x 


10" 


-3 


-1.1148 


x 


10" 


-4 


-4.3836 


x 


10" 




-5.2421 x 


10" 


-5 


4.9970 


P3 


3 


6.0592 


2.4699x10" 


-3 


1.9601 x 


10" 


-3 


-7.7253 


x 


10" 


•5 


-4.1596 


x 


10" 


-5 


-5.3403 x 


10" 


-5 


6.0635 


P4 


4 


7.0920 


2.1287x10" 


-3 


1.9532 x 


10" 


-3 


-5.7086 


x 


10" 




-3.5882 


x 


10" 




-5.1387 x 


10" 


-5 


7.0959 


P5 


5 


8.1662 


1.8559x10" 


-3 


1.8903 x 


10" 


■3 


-4.3311 


x 


10" 


■5 


-3.0363 


x 


10" 


■5 


-4.4437 x 


10" 


-5 


8.1698 


P6 


6 


9.2940 


1.6377x10" 


-3 


1.8491 x 


10" 


■3 


-3.3640 


x 


10" 


■5 


-2.6280 


x 


10" 


■5 


-3.8883 x 


10" 


-5 


9.2974 


P7 


7 


10.458 


1.4602x10" 


-3 


1.6973 x 


10" 


■3 


-2.6688 


x 


10" 


-5 


-2.1554 


x 


10" 


■5 


-3.4004 x 


10" 


-5 


10.461 


P8 


8 


11.641 


1.3166x10" 


-3 


1.1381 x 


10" 


■3 


-2.1654 


x 


10" 


■5 


-1.3052 


x 


10" 


■5 


-3.1208 x 


10" 


-5 


11.643 


P9 


9 


12.820 


1.2004x10" 


-3 


-4.7685 x 


10" 


-4 


-1.7971 


x 


10" 


■5 


4.9870 


x 


10" 


■6 


-3.0324 x 


10" 


-5 


12.821 


PlO 


10 


13.973 


1.1071x10" 


-3 


-4.3439 x 


10" 


■3 


-1.5262 


x 


10" 


■5 


4.1839 


x 


10" 


■5 


-3.1277 x 


10" 


-5 


13.970 


Pll 


11 


15.092 


1.0303x10" 


-3 


-1.1384 x 


10" 


-2 


-1.3205 


x 


10" 


■5 


1.0182 


x 


10" 


-4 


-3.2121 x 


10" 


-5 


15.082 


Pl2 


12 


16.197 


9.6297x10" 


-4 


-2.0660 x 


10" 


-2 


-1.1524 


x 


10" 


5 


1.7254 


x 


10" 


-4 


-3.0385 x 


10" 


-5 


16.177 


Pl3 


13 


17.312 


9.0177x10" 


-4 


-3.2143 x 


10" 


-2 


-1.0097 


x 


10" 


■5 


2.5156 


x 


10" 


-4 


-2.6884 x 


10" 


-5 


17.281 


Pl4 


14 


18.433 


8.4703x10" 


-4 


-5.1650 x 


10" 


-2 


-8.9005 


x 


10" 


■6 


3.8020 


x 


10" 


-4 


-2.3399 x 


10" 


-5 


18.383 



Table 5. Same as Table 2, tor f and p modes with (I, m) = (2, 2) and n = (0, . . . , 14). 



Mode 


n 


(TO 


,T T 
^2 




- 2 D 






,T T 

03 












-3° 




Ctot 


f 





2.7478 


6.7720x10" 


-4 


-1.4183 x 


10" 


3 


-5.3833 x 


10" 


■5 


8.8995 x 


10" 


-5 


1.0738x10" 


-4 


2.7472 


Pi 


1 


3.6591 


1.2040x10" 


-3 


-2.1590 x 


10" 


■3 


-8.7827 x 


10" 


■5 


1.3034 x 


10" 


-4 


1.7320x10" 


-4 


3.6584 


P2 


2 


4.8688 


1.1437x10" 


-3 


-2.2764 x 


10" 


-3 


-6.8308 x 


10" 


-5 


1.1444 x 


10" 


-4 


1.9652x10" 


-4 


4.8679 


P3 


3 


5.9318 


9.8595x10" 


-4 


-1.8483 x 


10" 


■3 


-4.9668 x 


10" 


-5 


7.8769 x 


10" 


-5 


2.2995x10" 


-4 


5.9312 


P4 


4 


6.9627 


8.5312x10" 


-4 


-1.2242 x 


10" 


■3 


-3.7113 x 


10" 


-5 


4.5164 x 


10" 


-5 


2.3593x10" 


-4 


6.9625 


P5 


5 


8.0358 


7.4678x10" 


-4 


-9.1572 x 


10" 


-4 


-2.8358 x 


10" 


-5 


2.9529 x 


10" 


-5 


2.1029x10" 


-4 


8.0359 


P6 


6 


9.1626 


6.6237x10" 


-4 


-6.1137 x 


10" 


-4 


-2.2221 x 


10" 


■5 


1.7440 x 


10" 


-5 


1.9074x10" 


-4 


9.1629 


P7 


7 


10.326 


5.9321x10" 


-4 


-1.0744 x 


10" 


-4 


-1.7761 x 


10" 


■5 


2.7379 x 


10" 


-6 


1.7293x10" 


-4 


10.326 


P8 


8 


11.508 


5.3505x10" 


-4 


1.3941 x 


10" 


■3 


-1.4449 x 


10" 


■5 


-3.2072 x 


10" 


-5 


1.6570x10" 


-4 


11.510 


P9 


9 


12.687 


4.8479x10" 


-4 


5.2270 x 


10" 


■3 


-1.1928 x 


10" 


■5 


-1.0965 x 


10" 


-4 


1.6851x10" 


-4 


12.693 


PlO 


10 


13.840 


4.4058x10" 


-4 


1.3779 x 


10" 


-2 


-9.9748 x 


10" 


■6 


-2.6615 x 


10" 


-4 
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Fig. 3. Frequency corrections against radial order n for p modes with (l,m) = (1,1) and 
n — (1, . . . , 15) for a zero-age main-sequence star with M = 12Mq. Here oo is the zero-order 
eigenfrequency, crj and aj are the second and the third-order Coriolis contributions, and erf 
and erf are the second- and the third-order non-spherically-symmetric distortions, a§ is the 
third-order distortion and Coriolis coupling, and a c = a J + aj + o% + erf + trf is the total 
third-order frequency correction. All frequencies are in units of \J GM/R 3 = 2.42 x 1CP 4 s _1 . 





Fig. 4. Same as Fig. 3, for p modes with (I, m) = (0, 0) and n = (1, . . . , 15). There is no aj , cr°, 
and a§ — for m = (see Eqs. (53), (65), and (69)). Also ctJ? vanishes since Z = m = and the 
model is a uniformly rotating star (see Eq. (64)). 




n (Radial Mode Number) 



Fig. 5. Same as Fig. 3, for f and p modes with (I, m) = (2, 1) and n = (0, . . . , 14). 



Zero Age Main Sequence Star (l,m)=(2,2) 




n (Radial Mode Number) 

Fig. 6. Same as Fig. 3, for f and p modes with (I, m) = (2, 2) and n = (0, . . . , 14). 



Table 6. Values ol zero-order eigenlrequency ao, total frequency a± corrected up to third order, 
total third-order frequency correction Aa± = (a± —ao) due to rotation and coupling, expansion 
coefficients A\ and A2 normalized to 1, i.e., A1+A2 = 1, for selected pairs of near-degenerate 
poloidal modes in a zero-age main-sequence star with M — 12Mq. All frequencies are in units 
of y/GM/R* = 2.42 x 1(T 4 s _1 . 



m 


/ 


n 


coupling 


CT() 




A<7± 




4 ±} 










1 


Pi 


2.9941 


3.0027 


8.5940x10" 


-3 


0.99950 


-0.03152 





2 





f 


2.9199 


2.9262 


6.2499x10" 


-3 


0.04412 


0.99903 








3 


P3 


5.1853 


5.1907 


5.3200x10" 


-3 


0.99719 


0.07496 





2 


2 


P2 


5.1189 


5.1266 


7.6955x10" 


-3 


0.08727 


-0.99618 








9 


P9 


11.932 


11.934 


2.5947x10" 


-3 


0.99858 


0.05319 





2 


8 


P8 


11.775 


11.782 


7.0084x10" 


-3 


0.05738 


-0.99835 



b. Concluding remarks 

The third-order effect of rotation on p and g modes for a uniformly rotating zero-age 
main-sequence star of mass 12 M Q has been investigated. The third-order perturbation 
formalism presented by Soufi et al. (1998) was used and revised because of some misprints 
and missing terms in some of their equations. Following Soufi et al. (1998), the Coriolis 
and spherically-symmetric distortion effects were included in the zero-order eigensystem. 
This yields eigenfrequencies ujo of eigenmodes which are no longer m-degenerate, even at 
zero order. Furthermore this procedure enables one to obtain eigenfrequencies with the 
required e 3 accuracy without the computation of eigenfunction corrections at successive 
orders of e. The zero-order eigenvalue problem was solved by pulsation code modified in 
this manner. Numerical calculations of oscillation frequencies were carried out for our 
selected model and second- and third-order frequency corrections due to Coriolis, non- 
spherically-symmetric distortion and Coriolis-distortion coupling were computed. For 
the case of m = there is no first and third-order frequency corrections. For the case of 
radial oscillation (I = m = 0) the second-order non-spherically-symmetric distortion is 
also zero and only the second-order Coriolis contribution exist. We discuss the validity of 
neglecting the surface terms which arise when the density derivatives are removed through 
an integration by parts. They become significant for higher-order modes, particularly in 
the present model whose atmosphere is relatively thin. Coupling only occurs between two 
poloidal modes with the same m and with I differing by or 2. 

We have carried out a careful comparison with the results of the independent imple- 
mentation of the third-order formalism by Soufi et al. (1998). After taking into account 
the modifications discussed in Appendix A the results of the two formalisms for the 
combined second- and third-order corrections agree to within a few per cent. 

In a subsequent paper we intend to investigate numerically the effect of rotation up 
to third order for a sequence of (i Cephei star models with uniform or radially varying 
rotation profiles. 
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Appendix A: summary ot modifications compared with the formulation ot 
S98 

The third-order perturbation formalism for the influence of rotation on stellar oscillations 
that is presented here is to a large extent a re-derivation of the results of S98. We 
consider it necessary to present some results of the derivations because of a number of 
typographical errors in the paper of S98, which might lead to confusion. For convenience 
we here summarize the required corrections to that paper, and other modifications in the 
present formulation: 

— There are several corrections to S98, Eqs. (A10). In the second of these equations the 
term — ftf/C is mistakenly omitted (cf. Eq. (21)). In the same equation U/(l — ay) 
must be used instead of U in the radial case (cf. Eq. (35)). In the last of S98, Eqs. 
(A10), the term in v should be 

UV g 



A- 



1 — <7 r 

(cf. Eq. (23)). 

— In the definition of (, S98 Eqs. (All), C r must be used instead of C (cf. Eq. (31)). 

— In their Eq. (A15), S98 neglected the factor (g c /g) as being 0(0 2 ) and approximated 
{g e /g)vk by Vk- We use the full expression in Eq. (41). 

— The factor 2J/uiq in Eq. (A22) of S98 was mistakenly omitted. The corrected form of 
the relation should be 

(£ti 9 | p£tifc) = (W.m, (%) i S h,i q K i+ ^ 

5l k ,l q + 2 K 2kq + Sl k J g - 2 K* qk ] , 

where K\ and K 2kq are given by Eq. (B5) in S98. 
A.l. Surface effects 

We have considered two approaches for the computation of the frequency corrections. In 
the first, the corrections are calculated from Eqs. (B6)-(B7), (B19)-(B21), and (B25) of 
S98. In the second which is used in present work, those coefficients are computed from 
Eqs. (52), (53), (64), (65), and (69). 

There are two substantial differences between the equations used in the two ap- 
proaches. In the formulation of S98, the density derivatives are eliminated through an 
integration by parts and the resulting surface terms are ignored. For instance, in Eq. 
(B25) of S98 the neglected surface term is 

T77 O 3 

S.T = -— r 5 u 2 (l + V )(S lklq9l + Q kq2 g 2 )p\ r=R , (A.2) 

where g\ — Ca^z^Zq and g 2 = C<T$(ykZ q + y q Zk + 3zkZ q ). In Fig. A.l, the surface term 
given by Eq. (A.2) as well as its effect on the third-order distortion and Coriolis coupling 



correction, o"3 Soufi (see Eq. {alb) m b98), arc plotted against radial order. Ihe plot 
shows that for lower radial orders n the effect of surface term is negligible whereas for 
higher orders, it becomes important. Note that for higher radial orders, the dominant 
term in surface term comes mostly from g 2 - Our study shows that for the model used in 
the present paper the surface terms become considerable. 

A. 2. Using the approximation z ~ y t /C<7Q in the third-order correction terms 

The other important difference which should be noted is that in the S98 approach for 
computing the third-order correction terms the approximation z ~ j/t/Cfg, which is 
valid for the non-rotating case, is used. In the present paper, on the other hand, the 
exact relation Eq. (24) is used. To investigate this difference in detail, in Fig. A. 2 the 
two terms included in Eq. (24) as well as the full relation for z are plotted for several 
p modes with (l,m) — (1,1) and n = (3,8,15) in the 12M Q ZAMS model considered 
here. Figure A. 2 shows that from the centre to near the surface the second term hiy/A 
is negligible compared with the first term C?/t/A. However, very close to the surface the 
second term becomes important and with increasing radial order n it becomes dominant 
compared with the first term. Hence we conclude that the approximation of neglecting 
the second term everywhere, particularly near the surface, is not valid. We find that the 
magnitude of this difference between the two approaches is more significant than the 
magnitude of the difference due to of the surface terms. 

If we include the surface terms in the S98 approach and use the approximation z ~ 
Ut/CdQ in the present formulation the results of the two numerical approaches are in good 
agreement. For lower radial orders, the difference in the combined second- and third-order 
frequency correction <j c is then smaller than one percent and increases slightly for higher 
orders. 





Fig. A.l. Surface term and its effect on Eq. (B25) in Soufi ct al. (1998) for p modes with 
(I, m) — (2, 2) and n = (0, . . . , 14) for a zero-age main-sequence model. 
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Fig. A. 2. The first Cj/t/A (left), the second hiy/A (middle), and the total terms (right) of zero- 
order horizontal component z of eigenfunction against fractional radius x = r/R for selected p 
modes (I, m) = (1, 1) and n = (3, 8, 15) for the model of Table 1. 
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Appendix b: Uitrerent contributions ot non-sphencally-symmetric distortion 

The dimensionless off-diagonal terms D kq , Eq. (62) are split into five separate integrals 
for convenience: 

D S = -^r- t d^ 4 riA fc A^ 2 , (B.l) 

Ukq JO 
r 2 /•! 



D( S = £r I d ^ 4 (! - °r){b 2 -(A + V g )u 2 } 

^Jkq JO 



kq JO 

X {Vq w k + WqHk + (A — i){z q V k + V q Z k )} , (B.2) 



U 2 



x[z q (A q -A k + 6)(X k -(A + V g )y k ) 

+ (\q-(A + V g )yq)(Ak-Aq + 6)zk} ( B - 3 ) 

+[b 2 -(A + V g )u 2 + Ul + m) + h% L ] 
x [Xky q + Vk\ ~2(A + V g )y k y q ] j , 

^"Jkq Jo 

x [6 2 - (A + V g )u 2 ] [(U + x- 4)y k y q 

+ (A g + A q z q )y k + (A fe + A k z k )y q ] 
-y k y q (A + V g )[(2 - a r ){b 2 -(A + V g )u 2 ) 

-2b 2 [y k y q + \y k z q (A q - A k + 6) 
+ il/,«fc(A fc -A, + 6)]} , 



x [6 2 - (A + V g )u 2 ] [y k y q + (A - 3)z fe z 9 ] , 



where A = (A fe + A q )/2 and 



Afe - 3m? 



(B.4) 



(B.5) 



(B.6) 



~ "4AfeT 

Appendix C: The hermiticity of an oscillating rotating system 

The hermiticity of the equations of oscillation in a rotating fluid has already been in- 
vestigated by Gough & Thompson (1990) who ensured the hermiticity for the rotating 
star by means of mapping each point in the distorted model to a corresponding point 
in the spherically symmetric stellar model. Lyndon-Bell & Ostriker (1967) demonstrated 
that all of the linearized operators appearing in a rotating system for a constant rotation 
profile are hermitian. Here we extend the argument of Lynden-Bcll & Ostriker (1967) to 
the case of a rotation profile that is a function of radius. 



lo study the hermiticity ot the operator L, Eq. (21) m b98, we need to show only 
that the operator A, Eq. (22) in S98, is hermitian. The other operators B, D and C 
which appear in Eq. (22) of S98, also appear in the total coupling coefficient H (see Eqs. 
(B1)-(B2) of S98), which is symmetric and hermitian, as is shown in Sect. 4. 

Let £ and ( be two displacement eigenvectors of an oscillating rotating system. Then 
one can show that the operator A, Eq. (22) of S98, satisfies the following relation 

<C|^) = <AC|0 + S.T., (C.l) 

where the surface term (S.T.) is 

s.t. = $ s {p'c + (C* • Vp)£ + (iw • C*)£ 

+p(4>'c - « (C.2) 

+ 4^(W*-^*V0')}-dS, 

with dS = e r R 2 sin 9A0d<p; here hat and tilde indicate components of the eigenfunctions 
associated with ( and £, respectively. In particular, <p' is given by the Poisson equation, 
Eq. (12), for eigenvector Q. Following Unno et al. (1989), the relevant boundary conditions 
for an oscillating rotating star are 

p = p = aXr = R, (C.3) 

p = at r = R , (C.4) 

^ + ^V = a tr = R. (C.5) 
dr r 

The first and second boundary conditions, Eqs. (C.3) and (C.4), result in removing all 
terms in Eq. (C.2) that include p and p', and also all terms that include Vp because 
they are proportional to p (see Eqs. (12) and (19) in S98). The last boundary condition 
eliminates the last surface term in Eq. (C.2). Hence all the surface terms go to zero at 
the surface of the star and therefore S.T. = 0, with the result that Eq. (C.l) becomes 

(C I AO = (A( | , (C.6) 

and therefore the operator A is hermitian. 

We note that from a theoretical point of view, the density vanishes exactly at the 
surface, whereas for a realistic model obtained from numerical calculations the density is 
not exactly zero at the surface. Hence from a numerical point of view the surface terms 
are not exactly zero; consequently, whether the surface term can be removed or not 
depends on the required accuracy. For instance, we saw before that the effect of surface 
terms is not negligible when the cubic frequency corrections are computed. 
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